JANUARY 1944 





of the 


UNIVERSITY OF PUERTO RICO 














~A Mernop ror THE So.ution or Norma Equations 
ut By Bernardo G. Capé 


Meruop or INTERPRETING THE REsULTs oF FIELD TRIALS 
By Bernardo G. Capé 


A New Mernop or PerrorMinG FIELD TRIALS 


4 By Bernardo G. Capé 


PUBLISHED BY 
THE AGRICULTURAL EXPERIMENT STATION 
Rfo Prepras, P. R. 





AGRICULTURAL EXPERIMENT STATION 
University of Puerto Rico Rio Piedras, Puerto Ri 


JAIME Benfrez, M.A., LL.M Chancellor 
T0168 OPERANT, WBS oo 65s aia Vice Chancellor 


STATION STAFF?! 


AnituRo Rous, MiBe soe. ko a a es Director 
R. Cotén Torres, M.S Asst. Director 
EB. DBL: TORO, IR. MBE. 6S oc oes ia 8h Admin. Asst. 
B. G. Capré, Pa.D Biometrician 
M. R. Garcfa Finance Officer 


Animal Husbandry: Phytopathology: ‘ 
, An. Husb., Head of Dept. L. A. Alvarez, M.S., Phytopath., Head of Dag 
José Adsuar, B.S., Assoc. Phytopath. 





Carlos M. Mufiiz, D.V.M., Assoc. Vet. 

F. J. Marchén, M.S.,? Asst. Chem. ae: 

Herman Cestero, B.S., Actg. Asst. An. Husb. Plont Breeding: P. 
, Plant Breeder, Head of De 


Victor L. Quifiones, M.S., Actg. Asst. Chem. OE oe ae 
A. Acosta Matienzo, M.S.,? Asst. Plant Breede 


Herbert R. G6mez, B.S.,* Actg. Asst, An. Husb. s 
Juan E. Veve, B.S., Actg. Asst. Plant Breeder! 


Agronomy and Horticulture: R. Vizcarrondo, B.S., Asst. Plant Breeder 

Arturo Riollano, M.S., Agron., Head of Dept. 
Pedro Richardson, B.S., Agron. Plant Physiology: 4 
Ferdinand Méndez, B.S., Assoc. Agron. in Charge ———_———, Plant Physiol., Head of D 

of Station Farm F, J. Ramfrez Silva, M.S., Assoc. Chem. | 
J. P. Rodriguez, B.S., Assoc. Agron. A. Gauthier Roig, B.S., Asst. Plant Physiol,” 
Policarpo Gonzélez, M.S8., Actg. Assoc. Agron. q 
J. Guiscafré Arrillaga, M.S.,2 Assoc. Agron. Soile: : 
Ce ee pater eee J. A. Bonnet, Ph.D., Soil Chem., Head of De 
L, A. Gémes, B.8., Asst. Agron. R. Rodriguez Torrent, B.8.,2 Assoc. Soil 
M. A, Diaz, B.S., Asst. Agron. A. Riera, B.S., Asst. Chem 
J. Lerfa Esmoris, B.S., Actg. Asst. Agron. P. Tirado Sdeeon BS. Jr. Soil Co: 

F. Mariota Trias, B.8., Jr. Soil Conserv. 


ae G. M. Soto, B.S., Jr. Chem 
Rafeel Arroyo, B.S., S.E., Chem., Head of Dept. » M. Soto, B.S., Jr. Chem. 


H. Cruz Monclova, M.S8., Ch.E.,? Assoc. Chem. Substations: 
J. H. Ramfrez, B.S., Assoc. Chem. 3 
Leonardo Igaravidez, B.S., Asst. Chem. Isabela: 
Fernando Marrero, B.S., 8.E., Asst, Chem. ——_—_—_——, Agron. in Charge 
z hy: C. J. Clavell, B.S.,2 Asst. Agron. 
een R. K. Giles, B.S., Actg. Asst. Agron. 
R. Colén Torres, M.S., Econ., Head of Dept. F. J. Julié, B.S.,2 Asst. Agron. 
L. M. Géigel, M.S., Assoc. Econ. R. A. Irizarry, B.S., 2 . i 
F, de Jess, M.S., Assoc. Econ. Oo Hane, ae ane 
8. Diaz Pacheco, B.B.A.,* Assoc. Econ. Plant Propagation Substation: - 
J. J. Serrallés, Jr., M.8., Assoc. Econ. J. 8. Simons, B.S., Asst. Agron. in Charge | 
M. Roses, M.S.,? Asst. Econ. 4 
M. Vélez, B.B.A., Asst. Econ. Animal Production Substation: , 
M. Herndndez, B.S., Actg. Asst. Econ. E, 8. Vilella, B.S., Asst. An. Husb. 
J. O. Morales, B.S.,3 Asst. Econ. 
M. Pifieroc, B.S., Actg. Asst. Econ. School of Tropical Medicine: 
J. D. Rivera Anaya, B.S., Asst. Econ. P. E. Serrano, B.S8., Jr. Chem. 
P. B. VA4zquez Calcerrada, B.S., Actg. Asst. Econ. J. Lopez Candal, B.S., Jr. Chem. 


hi 


Engineering: lilies: 
M. L. Vicente, C.E., M.E., Engineer J. Rivera Ruiz, B.L.S., Librari 
Entomology: 
G. N. Wolcott, Ph.D., Ent., Head of Dept. 1 As of date of issue. 
L. F. Martorell, Ph.D., Assoc. Ent. 2On military leave. 

F. Sefn, Jr., B.S., Assoc. Ent. 2 On leave. 





Entered as second-class matter September 14, 1934, at the post office at Rfo Piedras, Puerto Rico, u 
the Act of August 24, 1912 











; THE JOURNAL OF AGRICULTURE 
°OF THE UNIVERSITY OF PUERTO RICO 


sed peeetetly by the Agricultural Experiment Station of the University of Puerto Rico, for the publi- 
on 0! 





articles by members of its personnel, or others, dealing with any of the more technical aspects of 
scientific agriculture in Puerto Rico or the Caribbean Area. 








ol. XXVIII January, 1944 No. 1 


METHOD FOR THE SOLUTION OF NORMAL EQUATIONS 
By BERNARDO G. CAPO 


. Biometrician 


; INTRODUCTION 


‘The methods employed for the solution of a set of normal equations have 
gen modified in recent years to employ modern calculating machines to 
we fullest extent. Remarkable among these modifications is the abbre- 
ated Doolittle solution due to Dwyer (2), which requires a considerably 
Meduced number of entries in the table of solution of such a system of equa- 
i Dwyer (2, p. 443) claims that this method “‘is unquestionably the 
fethod which should be used by those seeking numerical solutions to the 
lassical problems of multiple and partial correlation.” Hotelling’s ex- 
rience with this and other methods has led him to state (4, p. 16) that 
fo the mass of least-square and other problems in which the inverse of a 
rix is needed, the best procedure appears to begin with one of the 
ne hods described by Dwyer,” etc. It is the purpose of this article to call 
tention to a method which is simpler, at least for relatively small'numbers 
f constants to be fitted, than Dwyer’s abbreviation of the Doolittle 
nethod. 
a THEORY OF THE METHOD 
» The method is based on Laplace’s development of a determinant, stated 
a Bocher (1, p. 26) as follows: 
’*Pick out any m rows (or columns) from a determinant D, and form all 
 m-rowed determinants from this matrix. The sum of the products of 
ach of these minors by its algebraic complement is the value of D.” 
) The application of this principle to the solution of a set of normal equa- 
ons will be explained now. 
Tet 
4 qQ bi Ci *** Mm 
Qa be C2 +++ Me 
soe - |, and let, fori <j < k, and with reference 
to Ajgg...n ’ 








* Ma 
aj b; S& Y — . 
Aij = ne = ab; — aj, 
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for example, Ags = . :| = dads — ashe ; 
5 
bi 
Bis =|.) * | = dees — Dyes 
ij b; C; ej it» 
for example, Bu = ni | = bya — bate ; 
4 
a; b; Ci a, bi Cy 
Ain = }a; b; c;|; forexample, Aims = |a¢ be C4]; 
a be C% as be Ce 
a; C 
Awj =| * “| = acy — af; 
oj a; ¢; ij iy 
for example, Ans = st = = Als — Agle ; etc. 





Now, the application of the above-mentioned principle will show that 
Ays = ACs — AisCe + Ars 
Ass = Atzsds — Ards + Aissde — Aouds , ete. 

The sign of any given term will be positive if “‘0”’ or an even number of 
inversions of adjacent subscripts is required to bring them into an ascend- 
ing order, and negative if said number of inversions is odd. 

Fisher’s (3, p. 144) modification of Gauss’ (5, p. 60 and 84) method of 
undetermined multipliers applied to the set of normal equations: 

Aa, + Bb; + sere ot Mm = Say, 
Ag + Bh, + --- + Mm = Sby, 


and Aa, + Bb, + --- + Mm, = Smy, 

where a, = Sa’, b; = a2 = Sab, --+ , bo = Sb’, c, = bs = She; ete., requires 
the solution of the following sets of equations for Caa, Cas, +++» Cw, Cre, 
eee, ond Cun: 


Ist set: Casdi + Cab) +--+ + Camm = 1 
Caste + Casbe + i + Camm = 0 


CacOn + Cada + +++ + Camm, = 0 
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METHOD FOR SOLUTION OF NORMAL EQUATIONS 


2nd set: Cadi + Chi + --- + Comm = 0 
Corde + Crabe +--+ Comme — 


| 
_ 


CavOn + Cdn + ae + ComMn = 0 


nth set: Cama: + Combi t --+ + Cnnm = 
Cum + Combo + hia + CmmMe = 0 


oO 


Coil + ComDn + MES + CnmMn Ne 


If one now lets R = 1/Ayzs...n, one can obtain the following relations 
for the c;;’s. 
From the Ist set of equations, one obtains 


C.. = Bog...nk, 
Ca = —Acoss...nk? ’ 
Cae = Aoso4...nlt, 


—_ = + Aogs...nk?; 
from the 2nd set, one obtains similarly, 


Cy = Ajos4...nle 
Coe = —Arsos...nk, 


Com = shy34...nk} 
and finally, from the nth set, 
Sans = Ajg3...(n—-yR. 


The values of the c;;’s may then be found by obtaining the products of a 
number of determinants of the order (n — 1) by the reciprocal of A13...n, 
which itself may be evaluated from the relation 


Ayps...n i Ays...(n-yMan ae Aygs...(n—2)nM na + ae +Aoss...nMy ’ 


the last term being positive or negative according to the respective number 
of inversions of subscripts required as mentioned above. The determi- 
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nants which are coefficients of the m’s in the last expression are the co- 
efficients of R in the relations for Cimm, C(m—1ym,*** » Cam above. 


APPLICATION OF THE METHOD 


The application of the method will be described in connection with the 
solution of a set of 3 normal equations involving 3 predictors. 
Normal equations: 


Aa, + Bh, + Co, = Say, 
Aa, + Bb. + Cc. = Shy, 
Aas + Bb; a Ces 


Scy. 
From the previous relations, 

Caa = BosR = (bocz — bsc2)R, 

Cap = —Adosll = —(aeC3 — a3C2)R, 

Cac = Aoshk = (a2b3 — azbe)R, 

Cw = Arsh = (ayc3 — a3¢,)R, 

Coe = —ApR = — (aybs — azb,)R, 

Coc = ApR = (aybe — arb)R, 

Az = Ajol3 — Aigo + Aogty , 


and R = 1/Aj3. 

To evaluate the c;,’s, therefore, one needs to know Ay, Ai3, A103, Ass, 
Aoos , and B2; to start with. The values of these expressions, which consist 
of differences between products of two terms each, as pg — 1s, may be ob- 
tained from the machine without recording any intermediate figures. 
From these values and the original data, A123 and its reciprocal, R, may be 
evaluated, and from these, by setting R in the machine, the values of the 
c;;s directly. The values of A, B and C and their standard errors may then 
be obtained in the usual way. 

The application of the procedure outlined above to the same data from 
Carver used by Dwyer (2) to illustrate his method follows. 


a, = 1.000 bi = 0.313 c, = 0.280 Say = 0.495 
a = 0.313 b, = 1.000 C2 = 0.652 Sby = 0.650 


a3 = 0.280 bs = 0.652 cs; = 1.000 Scy = 0.803 
Syy = 1.000 
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Av = dib — abi = 0.902031 C..= AnR = 1.758997 
Ajs = aybs — a3b; = 0.564360 Cre = —AisR = —1.100525 
Ans = C3 — a3, = 0.921600 Cy = Ajo3h = 1.797157 
Ass = dab — azbe = — 0.075924 Cx. = Aosk = —0.148055 
Aas = Acxl3 — Azle = 0.130440 Cs a= — Aogl? = —0.254363 
Bos = bocs — bse = 0.574896 Cow= Beyk = 1.121070 


Aps = Appl - Aj3Co a AogCy = 0.512810 
R = 1/Axs3 = 1.950040 
A = CaaSay + CaSby + CaeScy = 0.270706 
B = CuSay + CySby + Cr-Scy = 0.158521 
C = C,-Say + Cy-Sby + C.-Scy = 0.623846 


Reduced Syy = 1 — ASay — BSby — CScy = 0.262014. 


The information contained in rows 4 to 10 inclusive, of table VII, of 
Dwyer’s article (2), which required 37 entries (disregarding the zeros), is 
found in the last 12 rows of the preceding scheme, which contains but 18 
entries. In addition, the preceeding scheme contains the values of Ca, 
C», and Cy, which do not appear in said table VII, but are essential for 
the statistical evaluation of differences between A, B, and C. The calcu- 
lation of these values would require still 3 more entries in that table, 
raising the total number of entries to 40. 

For the case of 4 normal equations with 4 predictors the respective num- 
bers of entries would be 63 in the abbreviated Doolittle solution method and 
40 in the proposed method. Two additional advantages of the proposed 
method as against Dwyer’s abbreviated Doolittle solution are that the 
operations tend to be simpler and errors tend to be of less importance in 
the proposed method, since an error will not invalidate such a large pro- 
portion of the work performed subsequently to it as it does in the Doolittle 
technique when no check éolumns are included. 


SUMMARY 


A new method of solving normal equations is presented. A comparison 
is made of the application of the method to a multiple regression case used 
by Dwyer to illustrate his abbreviation of the Doolittle method. The 
proposed method for relatively small numbers of unknowns at least, is 
simpler than Dwyer’s method, considered to date to be one, if not the most 
efficient for the purpose. 
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A METHOD OF INTERPRETING THE RESULTS OF 
FIELD TRIALS 


By BERNARDO G. CAPO 


Biometrician 


In the comparison by means of field trials of a relatively small number of 
agricultural practices, the complete randomized block and Latin square 
designs are most commonly employed by the majority of research workers. 
When a relatively large number of such practices are to be tested, other 
designs are used, such as the incomplete randomized blocks, including the 
lattices, and designs which make use of check plots in one form or another. 
The reason for this attitude is that the effects of the differences in fertility 
between the different plots within any given block, file or column, as the 
case may be, are not eliminated from the estimate of the experimental error 
obtained by interpreting tests performed according to the above-mentioned 
designs by the simple methods of statistical analysis developed for use with 
said types of experiments. 

It occasionally happens, however, that field tests are performed in places 
where the fertility of the soil varies a great deal between spots relatively 
near to one another. Under these conditions, and especially with crops in 
connection with which relatively large plots must be used, the usual method 
of statistical analysis yields too high an estimate of error and it is impossible 
to determine any significant difference between the practices under trial. 
The possibility of the use of some alternate method of interpreting the re- 
sults of such experiments might conceivably be of benefit in a number of 
cases. It is the purpose of this article to suggest the use of a method 
applicable to cases of this kind which occur when relatively small numbers 
of practices have been tested. 

The method is based on the assumption of a different effect constant for 
every different pair of adjacent plots in the field when an even number of 
plots has been used. If an odd number of plots has been used, the previous 
assumption is made for all but 3 plots lying together in the field, for all three 
of which the same effect constant is assumed to hold. By fitting then a 
multiple regression equation to the results of such an experiment, it is 
possible to obtain an estimate of the effect constant of each of the practices 
tested. Mathematical statisticians certainly need no more information 
with respect to the way in which the method works than the one conveyed: 
in the three preceding sentences. The rank and file of field workers, how- 
ever, are not so fortunate in this respect. Due to this, the explanation 
of the method will be performed by following a numerical example through- 
out for a case where an even number of plots was used. The slight modi- 
7 
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fication needed for the cases of odd numbers of plots will be dealt with 


further on. 


Dr1acraM I 


Field arrangement of alfalfa experiment 
































50 40 30 20 10 
A H G ed E 
30 26 30 35 44 
49 39 29 19 9 
C J I H G 
37 36 29 30 34 
48 38 28 18 8 
E B A J I 
59 38 22 29 33 
47 37 27 17 7 

I F E D C 
60 42 35 24 30 
46 36 26 16 6 
G D C B A 
62 52 45 38 33 
45 35 25 15 5 
D A J I H 
57 33 39 40 48 
44 34 24 14 4 
H E D C B 
43 38 43 39 49 
43 33 23 13 3 
E I H G F 
39 36 41 45 57 
42 32 22 12 2 
F C B A J 
43 42 47 40 52 
41 31 21 11 1 

J G F E D 
48 46 50 52 56 





Diagram I shows the field arrangement of the plots of an alfalfa experi- 
ment performed at the Agricultural Experiment Substation farm at Isabela 


by Messrs. L. A. Serrano and C. J. Clavell. 





In each case, the first number 
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INTERPRETING RESULTS OF FIELD TRIALS 9 


is the plot identification number, the letter beneath it is the treatment 
identification letter, and the number beneath it is the observed yield in tons 
of green alfalfa roughage per acre obtained from 16 successive cuttings. 

In table I the treatments are described, and in addition, the total and 
mean yields obtained with each treatment are presented. 

Although in this experiment ten different treatments were tested, each 
one of the treatments may be expressed in terms of the amounts of nitrogen, 
phosphorus and potassium applied per acre. Thus, instead of trying to 
determine whether the observed difference in mean yields caused under the 
effects of any two given treatments is significant or not, one may try to 
determine whether the nitrogen applications have affected the alfalfa yields 


TABLE I 
Yields in tons of green alfalfa roughage per acre 

















| Treatments 
| Yields 

Letter | Units applied Lbs./A. } 
» 7. K NH: | P20s K:0 | Total | Mean 

| { | | 
A 0 | o 0 o | -o:|' @ | mm | ae 
B . i 3 3 Oo | 48 | 216 | 211 | 42.2 
C 0 1 4 0 | 48 | 288 | 193 | 38.6 
D ae 3 36 | 48 | 216 | 232 | 46.4 
E 2 1 3 72 | 48 | 216 | 228 | 45.6 
F a ee 2 36 «| «(48 |} 144 | 297 | 45.4 
G 1 1 4 36 | 48 | 288 217 | 43.4 
H —_— 3 2% | o | 26 188 | 37.6 
I ; ry 3 36 | (96 216 = «198 | 39.6 
J 2 2 4 72 96 | 288 204 | 40.8 





significantly or not; and similarly for the other two elements. In order to 
determine the effects of these elements, one may assume that the yield of 
any plot is the additive result of a series of terms as follows: 


Nibn oe Pidbp a Kb; + B; = 3%; (1) 


where Y; is the yield of plot 7; B; is the yield constant of the block of which 
plot cforms part; NV; , P; , and K; are the units of nitrogen, phosphorus and 
potassium applied to the crop in plot 7; and b, , b, , and }; are the respec- 
tive increases in yield caused by the application of each unit of nitrogen, 
phosphorus or potassium respectively. b,, bp, and b, are termed partial 
regression coefficients. The values of these coefficients are then the ones to 
be tested to determine whether the applications of each respective element 
have affected the yields significantly or not. 

Thus, if one assumes block 1 to consist of plots 1 and 2, and assumes 
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further that the yield constant of the plots in block 1 is B; , one may write 
the following equation for plot 1: 


b, + bp + 3b, + Bi = 56, (2) 


since plot 1 received treatment D, which consisted in the application of 1 
unit of nitrogen, 1 unit of phosphorus and 3 units of potassium. One may 
similarly write the following equation for plot 2: 


2b, + 2bp + 4b, + Bi = 52. (3) 


In a similar way, one may write down equations for the other plots, ob- 
taining a total of 50 equations, one for each plot. In these equations there 
would be the 28 unknowns b, , bp , bk , Bi, Bo, +--+ , Bos. It is impossible 
to determine values for these 28 unknowns which will fit the 50 equations. 
The next best solution, that of finding the most probable values of the con- 
stants, must be resorted to. Once these most probable values are deter- 
mined, one may estimate the yield of any plot by means of equation (1) 
above. The difference between the actually observed vield on any given 
plot and the value of said yield estimated by means of equation (1) is the 
error of estimate of the yield of that plot. Thus, for plot (1) the error of 
estimate would be: 


d, = 56 — b, — b, — 3b. — Bi, (4) 


since the actually observed yield was 56 and the estimated yield would be 
the value of b, + b, + 3b, + B, when the most probable values of these 
constants are substituted in equation (4). 

According to the principle of Least Squares, the most probable values of 
constants like the ones under discussion are those which render the sum of 
the squares of the errors of estimate a minimum. ‘Thus, if d, is the error 
of estimate of plot 1, dz is the error of estimate of plot 2, ete., those values 
of the constants b, , bp, b., Bi, B:, +--+, Bos which would make Sd’ = 
di + di + --- + dio a minimum, would be their most probable values, and 
therefore, the ones to be used in equation (1) for estimating the yield of any 
given plot. 


Now, 
di = (56 — b, — bp — 3b, — Bi)’. (5) 
Similarly, 
d; = (52 — 2b, — 2b, — 4b, — By)’. (6) 
In a similar way one may find the values of d3 , di, --+ , dio. 


In applying the criteria for a minimum to the expression Sd’ = di + 
di + --- + dio = (56 — b, — b, — 3b, — B,)* + (62 — 2, — 2, 
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4b, — Bi)? +--+ + (30 — Bas)’, one obtains 28 equations, each one of 
which is a partial derivative of Sd’ with respect to one of the constants 
b,, bp, be, Bi, Bo, +++ , Bos equated to zero. The equation obtained on 
finding the partial derivative of Sd’ with respect to B, and simplifying is 


56 — bn — bp — 3b, — Bi + 52 — 2b, — 2b, — 4b. — Bi = 0. (7) 
From this equation one gets 
B, = 54 — 3b,/2 — 3b,/2 — 7b,/2. (8) 


If one now substitutes this value of B,; in equation (4) above, and simplifies, 
one gets 


d; = 2+ b,/2 + b,/2 + bi/2 = (di — de) /2. (9) 

Substituting the value of B, in the corresponding equation for plot 2, ie., 
d, = 52 — 2b, — 2b, — 4b, — Bi, 

one gets d2 = —2 — b,/2 — b,/2 — b,y/2 = — (di — d2)/2. (10) 


It will be noticed that by means of these substitutions B, has been elimi- 
nated from the equations corresponding to the first two plots and, since 
these are the only two equations where B; occurs, from the whole set of 50 
equations. In a similar way all the other B’s may be eliminated from the 
other equations resulting finally in a system of 50 equations in but the three 
unknowns b, , b,, and b,. In this new system of equations, the right-hand 
side of the equation corresponding to plot 1 will be equal to the right-hand 
side of that corresponding to plot 2 except that the signs are reversed. 
Similarly the right-hand sides of the equations corresponding to plots 3 
and 4 will be equal except for the signs changed; and similarly for the re- 
maining pairs of equations. 

Now, the most probable values of the three unknowns b, , by , and b, are, 
as stated above, those which would make Sd’ a minimum. From equa- 
tions (9) and (10) above, it follows that d, = —d:, and therefore, dj = 
d; and —2d\d, = 2d; = 2d; = di + d3. Thus, (d, — d2)” = di — 2dyd, + 
d; = 2d; + 2d3. Asimilar relation may be found in the cases of the corre- 
sponding errors of estimate of the plots belonging to the other blocks. 
Therefore, 


(ch — ds)" + (ds — ds)” + +++ + (dig — do)” = 2d; + Qdz + ds + di 
+ +++ 4+ Qdig + 2dio = 28d”. (11) 


For simplicity of calculations Sd’ is calculated in the method under dis- 
cussion by calculating 2Sd° by the use of the left-hand side of equation 
(11), and dividing by 2. The numerical values of the coefficients of the 
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squares and products of the unknowns which appear in Sd’ may thus be 
conveniently found in the example under study by arranging the calcula- 


tions as they appear in table II. 


TaBLeE II 


Calculation of the numerical values of the coefficients of squares and products 
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Difference 

l dp |dx 
Plots 
a1 BS [1-2-1 
3-4 | F-B/ 1 Oj-1 
5-6 H-A | 1} 0} 3 
&7 | I-C 1] 1/—1 
10-9 | EG! 1 0j-1 
11-12 | E-A | 2} 1) 8 
13-14 | G-C | 1 0! 0 
15-16 | I-B | 1] 1] 0 
18-17 | J-D | 1) 1) 1 
20-19 | F-H | 0j 1/-1| 
21-22 | F-B | 1] 0-1 
24-23 | D-H/| 0} 1/ 0 
26-25 | C-J |-2/-1] 0 
27-28 | E-A | 2] 1] 3/1 
30-29 | G-I | Oj-1) 1 
31-32 | GC | 1] 0] O 
34-33 | E-I 1/-1] 0 
36-35 | D-A |! 1] 1] 3/1 
37-38 | F-B | 1] O/-1 
39-40 | J-H | 1| 2| 1/1 
41-42 | J-F 1/ 1] 2 
44-43 | H-B | 1-1] 0 
4645  G-D | 0] 0} 1 
47-48 | I-E |-1] 1] 0 
49-50 | C-A | 0} 1) 4) 
Sum.. 


Coefficients (= Sum/2)..... 
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(9) (10) 
dydx| dydy 
1; —4 
—1 8 
3 15 
—l 3 
—1 10 
6 24 
0 6 
0 2 
1 5 
0 0 
—1 3 
0 0 
0}; —12 
6 26 
0 0 
0 4 
0 2 
3 19 
—1 4 
1 10 
2 5 
0 4 
0 0 
0; -1 
0 0 
18 |133 
9 | 66.5 


(11) 
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(14) | (15) | (16) 
dr dxdy dy 
1 —4 16 
1 —8 64 
9 45 | 225 
1 —3 9 
1 {-—10 | 100 
9 36 | 144 
0 0 36 
0 0 4 
af 5 25 
1 —5 25 
1 —3 9 
0 0 4 
0 0 36 
9 39 | 169 
1 1 1 
0 0 16 
0 0 4 
9 57 | 361 
1 —4 16 
1 10 | 100 
+ 10 25 
0 0 16 
1 5 25 
0 0 1 
6 28 49 
7 | 199 {1480 
3.5] 99.5 | 740 





Columns (1), (2) and (6) have been filled from the information given in 


diagram I. 


In finding the differences corresponding to any given pair of 


adjacent plots, said differences have been taken in such an order as to yield 
Thus, since the yield of plot 2 was 
smaller than the yield of plot 1, the differences are to be found by subtract- 
ing the data of plot 2 from the corresponding data of plot 1. 


a positive value of dy in every case. 


cated by writing “1 — 2” in column (1). Since plot 1 received treatment D 
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and plot 2 received treatment J, the corresponding entry in column (2) is 
D-J. Theentry in column (6) is 4(= 56 — 52). The rest of columns (1), 
(2), and (6) is filled in a similar way. Thus, for plots 33 and 34, the order 
will be 34 — 33 since the yield of plot 34 exceeded the yield of plot 33. In 
columns (2) and (6) the corresponding entries are E — I and 2, respectively. 

After filling columns (1), (2), and (6); columns (3), (4), and (5) must be 
filled. In filling them use is made of table I. The entries corresponding 
to plots 1 and 2 in these columns are found as follows: dy = units of N of 
treatment D minus the units of N of treatment J = 1 — 2 = —1; and like- 
wise, dp = 1 — 2 = —1l, anddx = 3 —4 = —1. For plots 33 and34, 
where the corresponding difference between treatments is EH — I, dy = 
2-1l1=1dp=1-—2= —-l,anddx = 3—3=0. Inasimilar way, 
the rest of columns (3), (4), and (5) may be filled. 

Once columns (1) to (6) are filled, the corresponding entries in columns 
(7) to (16) may be made by performing in each case the operation indicated 
by the heading at the top of each column. Thus, for the entries corre- 
sponding to plots 1 and 2, dy = dydy = (—1)(—1) = 1; dydp = (—1)(—1) 
= 1;dwdx = (—1)(—1) = 1; dydy = (—1)(4) = —4; dp = dpdp = (—1) 
(-1) = 1; dpde = (—1)(—1) = 1; dpdy = (—1) 4 = —4;etc. Theen- 
tries corresponding to plots 33 and 34 are, likewise, as follows: dy = (1)(1) 
= 1;dydp = (1)(—1) = —1; dydx = (1)(0) = 0; dydy = (1)(2) = 2;ete. 

Once these entries are made, it remains but to add the entries of columns 
(6) to (16) and to divide each sum by 2. Since 


Sd = b7,Sn* + 2b,b,Snp + 2b,b,Snk — 2b,Sny + b;,Sp° 
+ 2b,b.Spk — 2b,Spy + bSk° — 2b.Sky + Sy*, (12) 
one can write the value of Sd’ for this case as follows: 
Sd’ = 14.50%, + 2(5.5)bubp + 2(9)bnby — 2(66.5)D, 
+ 10b;, + 2(8)byb, — 2(38.5)bp + 33.5b, — 2(99.5)b. + 740. (13) 


As previously stated, the most probable values of b,, bp, and bj are 
those which will make Sd? minimum. Among the mathematical require- 
ments for Sd” to be a minimum, the requirements stated by three following 
equations must be fulfilled: 


b, Sn? + b,Snp + b,Snk = Sny, (14) 
b,Snp + b,Sp’ + b.Spk = Spy, (15) 
and b,»Snk + bpSpk + b,Sk’ = Sky. (16) 


Equations (14), (15), and (16) are known as “normal equation in b, ,”’ 
“normal equation in b, ,’’ and “normal equation in b; ,” respectively. In 
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the sample under study, the equations are as follows: 


14.5b, + 5.5b, + 9b, = 66.5, (17) 
5.5b, + 10b, + 8b, = 38.5, (18) 
and 9b, + 8b» + 33.5d, = 99.5. (19) 


These three simultaneous linear equations in 3 unknowns may be solved 
by any of the well-known methods studied in the elementary courses in 
algebra to obtain the values of the partial regression coefficients b, , b,, 
and b.. Due, however, to the fact that the knowledge of the standard 
error of any given statistic is almost as essential as the value of the statistic 
itself for the proper evaluation of the significance of said statistic and of 
differences between it and some other statistic, Fisher’s modification of 
Gauss’ method of correlatives or indeterminate multipliers is usually re- 
sorted to in practice. Snedecor (4, p. 302) and Rider (3, p. 39) present 
discussions and give examples of solutions of these systems of equations 
by this method. Doolittle’s method of solution of normal equations as 
discussed by Mills (2, p. 656) may be incorporated to advantage to the 
method of indeterminate multipliers. The author has found these methods, 
however, too tedious and intricate for use in practice. He suggests, there- 
fore, the use of the following method which, for a small number of con- 
stants to be fitted, is much easier and takes considerably less time to apply 
than the methods previously referred to, specially if one has a printed or 
mimeographed form stating the calculations to be performed, so that one 
has to fill in merely the values asked for in said form. The suggested 
method of solution is as follows: 

Following Fisher, (1, p. 144), the equations to be solved are: 


Caathh + Cadi + CoC = 1, (20) 
Conte + Cabs + Cale = 0, (21) 
Coats + Cabs + CacC's = 0, (22) 
Cat: + Cur + Cody = 0, (23) 
Cate + Car + O05 = I, (24) 
Cad; + Cubs + CC 3 = 0, (25) 
Cacti + Coed + CecCi = 0, (26) 
Catia + Code + Cale = 0; (27) 
and Cacds + Crcdbs + C.cCs = 1, (28) 
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where a; = Sn’, a2 = b; = Snp, ag = 1 = Snk, bk = Sp’, bs = c = Spk, 
and cs; = Sk’. 


: ’ b 
Now, if Ay represents the determinant of the second order js a ; 
| 42 U2 
a by C1 | 
Ay; means the determinant of the third order | az bz c:|, and, in general, 
| Ag bs C3 | 


i 


Aws...n Ieans the determinant of the nth order 


ay by +: my, 
| dz bo +++ Me | 
| 
| | 
| ' 
: | 
| | 
| 


Anda *** My | 


the following relations may be easily deducted from the fact that any given 
determinant is equal to the sum of all terms which may be formed by writing 
down the principal diagonal a,boc; --- m, and forming all possible permuta- 
tions of the subscripts; the sign of each term being positive if an even num- 
ber of inversions must be performed to bring the subscripts into an as- 
ending series and negative if said number of inversions is odd. 


Ay = dibs — ash 2 Ajo3 = AyoC3 — Asse + Asst ) 
Ajoss = Ajozdy — Arosds + Assi, — Asad , ete. 


If now Ayos = QyC3 — Asi ; Acos = Ae2C3 — AgCo and Boz = becs — bsc2, the 
values of Cac, Cas, and Ca. may be determined from equations (20), (21), 
and (22) above to be as follows: 


Coa = BosAis; Cos = —AopAiee; and Cac = AvsAizs. 
From equations (23), (24), and (25), the values of Cy and C;, come out to 
be as follows: 
Cw = AwAns; and Cy = —AAns. 
Finally, from equations (26), (27), and (28), 
Cee = AvAis. 


The number of operations required to evaluate these statistics is not 
large and does not take too long, specially if, as already stated, one has some 
form indicating the operations and has to fill only the numbers in it. Such 
aform, already filled, appears as Table III. 

In this table, A corresponds to b, , B corresponds to b, , and C corresponds 
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to b, which are the unknowns whose values were to be determined. Thex 
values of b, , bp , and b, may be checked by substituting them in equatioy 
(17), (18), and (19) above. 


TaBLe III 
Solution of simultaneous linear equations 





Aq, + Bb; + Cc, = 
Adz + Bb» + Cc 
Aa; + Bb; + Coy = 


a, = 14.5 b, = 5.5 

ag = 5.5 bo = 10 

a3 = 9 bs = § 

a,b. = 145 dob, = 30.25 
ab; = 116 Azbs = 44 
a\c3 = 485.75 Qsc3 = 184.25 


Aj. = ayb2 — axbi = 114.75 


Ai = abs aad a3b1 = 66.50 
A103 = A)1C3 — A3C0, = 404.75 
Ag; = debs — asb2 = —46 
Azo3 = G2C3 — AsC2 = 112.25 
Bas oe boc — dace = 271 
Coa = BA, = 0.0935088 
Cas = —Ar04 |, = —0.0387320 
Cac = AnA. = —0.0158723 
123 
Checks: 
Goats = 1.35588 C opQ2 = 
Cavb; = —0.213026 Crvbe = 
Cacti = —0.142851 Crt: = 
1.0000 
CaaSay = 6.218335 CawSay = 
CaSby = —1.491182 CySby = 
CaScy = —1.579294 CrScey = 
A = 3.147859 B= 


Say 
Sby 
Scy 


= —532 


Say = 66.5 
Sby = 38.5 
Scy = 99.5 
bees = 335 
b3Ce = 64 


cy = 9 

C2 = 8 

= 33.5 

a3b, = 49.5 

a3b2 = 90 

asc; = $1 

a3C2 = 72 

Aioc3 = 3844.125 
—Aj3C2 = 

Ao3C1 = —414 


Aj93 = 2898. 125 


12 


Cop - 
Coe - 


. 213026 
1.39659 
. 183567 


1.0000 
.575678 
5.376872 
2 


0.518077 


eat 
—A;:A,, 


A ; = 0.000345051 


3 


C. cz 
Cy cbs - 
Ce ccs = 


CacSay 
CrSby 
CeScy 


C 


A147, = 0.139659 
= —0.0229459 
= Andy, = 0.0395946 


= —0.142851 


—0.183567 


1.32642 


1.0000 


—1.055508 
—0.883417 


3.939663 


2.000738 


The significance of these coefficients of regression may be checked by 


calculating their standard errors and using the “t-test.” 


The total sum of 


squared deviations corrected for variations in fertility from block to block 
was found in table II to be 740, subject to 25 degrees of freedom, since 
from the total original number of 49 degrees of freedom, 24 (= 25 — 1) 


belong to the block statistics assumed to exist. 


The reduction in this 
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sum of squares caused by fitting the statistics b, ; b, , and by is found, as 
usual, by means of the relation 


Sy” = b,Sn” + b,Sp” + b,Sk’ = 209.33 + 19.95 + 199.07 = 428.35, (29) 


corresponding to 3 degrees of freedom, since it corresponds to the three 
statistics so fitted. The reduced sum of squared deviations is, therefore, 
s@ = S(Y — Y’) = 740 — 428.35 = 311.65, corresponding to 22df. 

The estimate of the variance, V, is, therefore, V = 311.65/22 = 14.1659. 
The variance of b, is then CaaV = (0.0935088)14.1659 = 1.3246. The 
standard error of b, = S.E., = (1.3246) = 1.151, and the corresponding 
value of ¢ = 3.1479/1.151 = 2.73, which is significant at the 5% point; 
since the value of ¢ at the 5% for 22df. is 2.074. 


TaBLe IV 


Analysis of the total sum of squared deviations 




















| ‘i | a ane F values 
5 tie deviati egrees | A Jariance 
ource OF the Ceviahions of freedom | Po ao estimate Experi 5% 1% 
eT | 49 | 4,651 
ae | 4 | 902 
Treatments............. a 922 | 102.4 | 1.26 2.15 2.94 
MR Sin senile 3 cuits | 86 | 2,927 | 81.31 | 











Similarly, Vz,.= 1.9784, S.E., = 1.407, and &, = 0.37, which is not 
significant. For b; , the corresponding figures are: Vo, = 0.5609, S.E., = 
0.7489, and &, = 2.67, which is also significant at the 5% point. 

The performed analysis indicates that the crop responded significantly 
to the applications of nitrogen and potash, whereas it did not do so with 
respect to those of phosphorus. 

The usual analysis of variance of the yield data of the experiment under 
study yielded table IV, the first block consisting of plots 1 to 10, the second 
block of plots 11 to 20, ete. 

The usual hypothesis thus fails to show any significant effect of the 
treatments on the yields. However, since in the previous analysis but 3 
statistics were fitted to explain the effects of the fertilizer applications 
whereas in this analysis 9 such statistics were fitted, a new analysis was 
made in an attempt to explain the sum of squares due to treatments in this 
last analysis by the use of only 3 statistics. Such an attempt indicated 
that of the sum of squares due to treatments, 922, a total of 334 could be 
explained by fitting the 3 statistics, there remaining 588 to be assigned to 
interactions of one sort or another. If this remainder is pooled with the 
sum of squares due to error, 2,927, a total of 3,515 is obtained, subject to 
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36 + 6 = 42df. The new estimate of the error variance would be, there- 
fore, 83,69, of the same order of magnitude as the value formerly obtained 
of 81.31. On testing the fitted statistics, whose values came out to be: 
b, = 2.167, b, = 0.529 and b, = 1.188, it is found that none of them js 
significant. 

On comparing the variance estimate obtained by the last method of 
analysis, 83.69, with that obtained by the proposed method, 14.1659, one 
obtains a relative precision of 5.9078, (= 83.69/14.1659), in favor of the 
2-plot blocks hypothesis as against the usual hypothesis, in this case the 
10-plot blocks hypothesis. There has been, in this case, therefore, an 
increase in precision of 491 per cent due to the change in hypothesis. — 

The results of this experiment might have been also interpreted by as- 
suming 5-plot incomplete blocks consisting of either the plots receiving 
treatments A, C, E, Gand J or treatments B, D, F, H, and J. The results 
of such an analysis appear in table V. 


TABLE V 


Analysis of the total sum of squared deviations 

















| Sunk ok F values 
Source of the deviations na. squared ; pa ————- 
deviations Experi- | 5% 19% 
mental | 0 ld 

Meee eee | 49 | 4,651 | 
BEUNONS 56 5i0156:aie «5 c.4ie-s'si as | 9 | 2,817 | 
Treatments.............] 8 830 103.75 2.21 | 2:25 3.12 
9s a ae ee | 32 1,504 47.00 | |" 








The observed differences between the treatment means are again not 
significant, although there has been an increase in precision of 81.31/47 — 
100% = 73%, by the use of the smaller 5-plot blocks. 

In case that an experiment consists of an odd number of plots, either by 
design or accident, the contribution of all but 3 plots to the numerical 
coefficients of the squares and products which appear in Sd” will be cal- 
culated as described above in the case of table II, by finding the respective 
differences between the adjacent plots, finding all possible squares and 
products of these differences, adding these squares and products, and 
dividing each of these sums by 2. 

The contribution of the 3 plots left out in the previous calculations 
may be computed as follows: The set of differences corresponding to the 
first of these plots will be formed by subtracting the sum of the corre- 
sponding figures of the second and third plot from twice the corresponding 
figures of the first plot. The set of differences corresponding to the second 
plot will be likewise formed by subtracting the sum of the corresponding 
figures of the first and third plot from twice the corresponding figures of the 
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second plot. In a similar way are found the differences corresponding 
to the third plot. 

Again the sums of all possible squares and products of these sets of 
differences are found, but in this case the sums must be divided by 9 in 
order to determine the contributions of these last 3 plots to the numerical 
coefficients. The sum of the contributions of both sets of plots will give 
the required numerical coefficients of the squares and products which ap- 
pear in Sd”. 

To illustrate the way of finding the contribution of these last 3 plots in 
such an analysis, one may assume that plots 48, 49, and 50 of diagram I 
are the remaining 3 plots in the case of an experiment with an uneven num- 
ber of plots. The value of dy for plot 48 would be found as follows: Since 
plot 48 received treatment L, with 2 units of nitrogen; plot 49 received 
treatment C, with no nitrogen; and plot 50 received treatment A, with 
no nitrogen, the respective difference, dy, would be, therefore, 
2/2) -0 — 0 = 4. The corresponding dp would be 2(1) — 1 — 0 = 1; 
the corresponding dx would be 2(3) — 4 — 0 = 2; and the corresponding 
dy would be 2(59) — 37 — 30 = 118 — 67 = 51. Ima similar way, the 
differences corresponding to plots 49 and 50 are found. 

In the above discussion, the question of the significance of the difference 
between any pair of statistics fitted to data of the kind mentioned has 
been dispensed with. In the cases where several varieties or field practices 
have been tested, the interest centers, not on whether a given variety or 
practice produces a significant departure of the measured effect from the 
mean effect of all the varieties or practices tested, but on whether the given 
variety or practice causes a significantly greater or smaller effect than some 
other variety or practice included in the test. In such cases, the number 
of statistics fitted to the data is one less than the number of treatments 
under study. Thus, if, say, varieties A, B, C, D, and E were tested in one 
of such experiments, statistics for A, B, C, and D only would be fitted, it 
being understood that the statistic for E is equal to the negative sum of 
the statistics corresponding to the other four varieties, ie, E = 
-(A+B+C+D). The corresponding C.,’s would be found by means 
of the relations: 


Cea = —(Caa + Cas + Cac + Caa), (30) 
Ca = —(Ca + Cru + Cre + Cra); (31) 
Cee = —(Cac + Cre + Cee + Cea), (32) 
Cea = —(Caa + Coa + Cea + Caa), (33) 
and Coe = —(Cac + Cre + Cee + Cie). (34) 


The variance estimate of a difference between any two such statistics, 
say A and B, would be found by means of the relation (Caa — 2Cas + Crs) V, 
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from which the corresponding values of the standard error and ¢ would be 
found accordingly. 

It now remains to present a resumé of the results obtained to date in the 
use of the 2-plot blocks method in the interpretation of the results of field 
trials. It may be stated that, in general, and in full accord with the 
theoretical foundations of the method, it is most useful as compared with 


TaBLeE VI 
Comparison of efficiency of the 2-plot blocks method with the usual analysis of 
“variance”? method of interpreting the results of complete randomized blocks 
experiments 























oe Increase in 
| variability | {ficiency 

Crop Nature of test | Factor studied ee by! use of 

a the 2-plot 

| analysis of | ios 

| method method 

Corn, Ist crop.......... Fertilizer | Yield 11.04 13.15 
Corn, 2nd crop. ...... 5 vs es 17:13 —6.43 
Com, 8rdcrop.......... - eS 9.97 | —11.00 
Alfalfa, 16 cuttings..... ‘6 unlimed | ‘“ 21.93 490.81 
Jd ts | oar €¢ limed | ‘“ 713 32.60 
MOOMNIO GER oos)5 ox 2 oa ake ne fas 16.79 27.86 
Sweet potato, lst crop.. ¢ 6 43.71 248.75 
4 jd 7°: Vad cae - S: 26.92 51.29 
Cueumbers, 1932-33... .. ee cs 18.38 11.31 
Cucumbers, 1933-34..... es Ke 18.20 109.32 
“ , 1935-36... .. “ “ 27.67 | 74.88 

“ , 1936-37... «“ “ 91.02 | 75.15 

«“ , 1941-42... «“ « 25.25 24.81 
OSES SE a ae is f 16.27 —7,61 
Sugar cane, Aguirre..... Varietal Yield of cane 17.27 53.48 

1942-43 
ms bale oe es Tons sugar/A 17.31 56.93 
ss e.7# La “ Sugar % cane 4.62 | —13.18 
Sugar cane, Rio Piedras 

“CULE ES Ug 2 ae ape J. Tons sugar/A | 27.60 190.16 





the usual method of analysis in cases where soil heterogeneity affects mark- 
edly the effect under study, the more so the more variable the soil. In 
table VI, the results of comparisons of the 2-plot blocks method with the 
usual method of analysis are presented. In every case the comparisons 
between the two methods have been made by using the same number of 
statistics to explain the effects of the treatments. 

As seen above, in by far the majority of the cases, the .2-plot blocks 
method has been more precise than the complete randomized blocks 
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method. The method, however, takes more time and work to apply to 
the results of trials not specially designed for its application. By using 
some sort of balanced design, however, the work of calculation of the 
results of experiments by means of the proposed method may be consider- 
ably shortened. 

Thus, in testing 6 treatments, each of the treatments might be included 
with each of the other ones an equal number of times in the 2-plot blocks. 
This would require replicating each treatment either 5 times or some 
multiple of 5. In general, this would require m(n — 1) replications where 
nis the number of treatments tested and m is any positive integer. Sucha 
layout for 6 treatments replicated 5 times could be something as illustrated 
in diagram IT. 


DiacraM II 


Layout for a plot distribution to test 6 treatments (A, B, C, D, E, and F) each 
replicated & times 
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The advantage of such a balanced design would be that all C;,’s would be 
equal, i.e., Coa = Co = Coo = Caa = Coe = Cy, and, similarly, all C;;’s 
would be equal, i.e., Ca = Cac = --- = Cez. Therefore, the same variance 
estimate would be used to test each effect constant and, likewise, there 
would be but one variance estimate to test the difference between any two 
such effect constants. More information relative to the possibilities of the 
use of these balanced designs will be given in another article to be published 
shortly after this one in this same Journal. In said article, the manuscript 
of which has been already prepared, the modified calculational technique 
is discussed and a numerical example is presented. . 
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A NEW METHOD OF PERFORMING FIELD TRIALS 


By BERNARDO G. CAPO 


Biometrician 
INTRODUCTION 


In a recent article (1) in which the author discussed a new method of 
interpreting the results of field trials, he presented evidence tending to 
show that the assumption of a different effect constant for every pair of 
adjacent plots of a randomized block field trial led to a greater precision 
in the reduction of the corresponding data than the assumption of a dif- 
ferent effect constant for every different complete block of the experiment. 
In said article he presented a method which might be used in this connec- 
tion in order to reduce the calculational work to a minimum for the inter- 
pretation of such experiments by use of the proposed new assumption. 
The work involved in the application of said method to the results of 
randomized block experiments already performed, however, is greater than 
the work required for the application of the usual methods of statistical 
analysis in regard to which said experiments have been designed. On the 
other hand, the increase in precision obtained by the application of the 
suggested method seems to the author to warrant fully the increase in such 
calculational work. 

As it was also mentioned in the article referred to, it is possible to reduce 
the number of operations required for the application of that method to the 
interpretation of field trials by designing such field trials so as to make use 
of the proposed method. The purpose of this article is to suggest designs 
for possible use in this connection as well as to discuss the procedure to be 
followed in the interpretation of such specially designed experiments. 


REQUIREMENTS TO BE FULFILLED BY FIELD TRIALS 


Field trials must fulfill certain requirements in order that the application 
of the proposed method to the interpretation of their results be most eco- 
nomical in time and effort. These requirements are as follows: 1—As many 
different groups of two treatments each must be made as is possible. 2— 
The experimental field must be divided into blocks of 2-adjacent plots each, 
so that the topographies of the two plots in any given block will be as nearly 
alike as possible. 3—Each group of two treatments must be assigned at 
random to be tested in one of the 2-plot blocks. The plot of the respective 
2-plot block in which one or the other of the two treatments of the group 
assigned to it will be tested must be also decided at random. 4—The same 
number of replications of each group of treatments must be used. Thatis, 
22 
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if n treatments are to be tried, there should be m(n — 1) replications, 
where m is any positive integer. Thus, if 6 treatments are to be compared, 
either 5 or 10 or 15, ete. replications of the treatments are required. 


EXAMPLES OF DESIGNS OF POSSIBLE USE 


In order to give illustrations of the possibilities of this type of layout, 
two examples will be given of possible different designs for a trial in which 
5 treatments are repeated 4 times each. 

Where 5 treatments, say A, B, C, D, and E, are to be repeated 4 times 
each, the experiment should include 10 groupings as follows: AB, AC, 
AD, AE, BC, BD, BE, CD, CE, and DE. The 20-plot experimental 
field is divided into ten 2-plot blocks, and the above treatment pairs are 


Di1acraM I 
Possible layouts of plots 
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DB.CD 
assigned at random to be tested in one or another of the 10 blocks, there 
being also randomization of the order of the treatments to be tested within 
each block. A possible geographical distribution might be as indicated 
by Diagram I (a). 

Another possible arrangement, which lends itself to possible interpre- 
tation by either the 2-plot block method or by the ordinary method of 
analysis for a 4-randomized-block trial, is shown in Diagram I (b). 


THEORY OF THE PROCEDURE OF CALCULATION 

As indicated in the article (1) previously referred to, the regression 
equation to be used in explaining the results of such an experiment would be 
Y,;= B+ T;; (1) 

where Y; is the observed effect, say yield, of plot 7, B; is the effect constant 
of the block of which plot 7 forms part, and 7’; is the effect constant of the 
treatment tested in plot 7. After writing down the respective equation 
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for each of the plots used, the problem becomes one of finding the most 
probable values of the B,’s and the 7’,’s, and the standard errors of the latter 
for comparison with one another. It was suggested in said article that a 
convenient way to do this would be one in which the differences between the 
corresponding figures of the paired plots were found and then, following 
the usual methods of regression analysis, determining the most probable 
values of the 7,’s and their standard errors. 

Now, if the experiment has been carried out in such a way that every 
treatment has been thus paired with every other treatment once, it will be 
found that by the application of the above procedure, the value of each of 
the C;’s will come out to be 2(n — 1)/n’ and the value of each of the C;/s 
(i ¥ j) will come out to be —2/n’ where n is the number of treatments 
tested. 

By knowing then Say, Sby, Scy, Sdy, Sey, and Syy, and by the proper 
use of the C;,’s and C;;’s, already known, the values of the treatment effect 
constants, i.e., A, B, C, D, and FE, may be found together with their stand- 
ard errors. One may thus test the significance between the estimated 
difference in effects between any two of the treatments tested. 

The theory of the procedure outlined in the last paragraph is so well 
known to those interested in it, that it does not merit further consideration. 
It may be worthwhile, however, to derive the relations C;; = 2(n — 1)/n 
and C;;(i ~ j) = —2/n*. This will be done now. 

Let it be assumed that n treatments have been tested in a trial of the 
kind under consideration. Let the effects of the different treatments above 
or below the mean effect of all treatments be designated by A, B, C, --- , N, 
where JN is the effect of nth treatment. The effect equations of the type 


Y;=B;+T; 


may then be written for each plot. 

Since any given treatment, say A, has been replicated n — 1 times, there 
will be n — 1 such equations in which the treatment effect A appears. Also, 
on account of this, there will be » — 1 equations formed by finding the 
differences between the equations corresponding to the treatment A plots 
and those corresponding to the other plots paired to said treatment A 
plots. These equations might be written as: 


+A + B = dsaze 


+A + C = dsaxc 
+A + D = ds.azp 


+A + (N — 1) = dsazw-» 
+A + N = dian 
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where the d’s stand for the respective differences between the effects noticed 
in the A plots as against the effects noticed in the plots paired with the A 
plots. Since, however, A, B, C, --- , and N have been assumed to be the 
effects of the treatments above or below the mean effect of all the treat- 
ments, 


A+B+C+-:-+N=0, 
or 
N=-(A+B+C+---+(N —1)). 


The last of the above equations, therefore, should appear rather as 
+2A+BaCxa+-:-:: + (N — 1) = dyazgy, due to the relation for N 
previously pointed out. Further, the effect constant A would also appear 
in all the equations corresponding to comparisons between plots receiving 
treatment N and anyone of the other treatments. 

In other words, A appears with a coefficient of +1 in the n — 2 equations 
corresponding to treatment A comparisons with treatments B, C, D,--- , 
N — 1, and also in the n — 2 equations corresponding to the comparisons 
of treatments B, C, D, --- , and N — 1 with treatment N; and furthermore, 
with a coefficient of +2, in the equation corresponding to its comparison 
with treatment N. That is, A will appear with a coefficient of +1 in 
2n — 4 equations and with a coefficient of +2 in one equation. This is 
true also of each of the other treatment effect constants. 

Therefore, the sum of the squares of the A _ coefficients will be 
(2n — 4) + 4 = 2n. This is true also of the sum of the squares of the 
coefficients of any other treatment effect constant. But since, as was 
demonstrated in the former article, on finding the sums of squares and 
products of the differences between the coefficients of the treatment effect 
constants in the original equations, one would be finding really 2Sa’, 
2Sab, --- , 2Sa(n — 1) instead of Sa’, Sab, --- , Sa(n — 1) which are 
the figures sought, the above value of 2n must be divided by 2, so that the 
value of Sa? = Sb’ = --- = S(n — 1)? =n. 

Now, in the equation corresponding to the comparison of A with B, one 
of the treatment effect constants is +1 and the other —1, so that their 
product is —1. In the equation corresponding to the comparison of A with 
N, the coefficient of A will be +2 and that of B will be +1, so that the 
products of these coefficients will be +2. This is true also of the equation 
corresponding to the comparison of B with N, where the coefficient of B 
is +2 and the coefficient of A is +1. In each of the n — 3 equations 
corresponding to the comparisons of each of the other treatments with N, 
the coefficients of both A and B will be +1, so that their product in each 
case will be +1. Therefore, the sum of the products of the coefficients of 
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A and B in these equations will be —1 + 2+ 2+ (n —3) =n. Again, 

this must be divided by 2 in order to get the proper value of Sab. Thus, 

Sab = Sac = --- = Sa(n — 1) = She = --- = S(n — 2)(n — 1) = n/2, 
The normal equations will be, therefore, 


An + Bn/2 + Cn/2 + --- + (N — 1)n/2 
An/2 + Bn + Cn/2 + --- + (N — 1)n/2 


Say, 
Sby, 


An/2 + Bn/2 + Cn/2 + --» + (N — 1)n = S(n — l)y, 


giving rise to the (n — 1) columned determinant 








n n/2 n/2 +--+ n/2| 
n/2 n n/2 +++ n/2 
D, = 
n/2 n/2 n/2 --» ni 
(n — 1) columns 
from which the values of Caz, Cas, +++, Con—1y(m—1) May be evaluated as 


follows. 

By subtracting the last row in D, from each of the other rows, adding 
the first to the last column in the expression thus obtained, expanding the 
expression thus obtained by minors along the first row, and by continuing 
such additions of the first to the last column and expansions by minors 
along the first row in the resulting expressions, one obtains finally 


- 0 
n n-319 n n—2 n = 2\" 
p, = (3) n oar - (3) [n+ m2) 2] =n(3) 5 
2 2 


Now, Caa = Naa/D, , where 
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n/2 n <-++ n/2 
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|n/2 n/2 +++ Mm | 





- 
(n — 2) columns 








Goi 
press 
No-t 





ain, 
lus, 
1/2, 


ing 
the 
ing 
ors 


1 











NEW METHOD OF PERFORMING FIELD TRIALS 27 


The evaluation of Naa, following a procedure identical to that used in 
evaluating D, , yields the expression 


na | n/2 0 
Noo = ) n 
— (5 n/2 n+ (n— 3) 5 


= ey E + (n — 3) =| = (n — 1) (3). 


(‘xa is, therefore, equal to 2(n — 1)/n”. On finding the corresponding ex- 
pressions for the other treatments, it is found that Ni = Ny =--- = 
Nn—1)(n—1) and, therefore, Cw ad Cw eee eS C(n—1) (2-1) = 2(n —_— 1)/n’. 


Now, Ca = Na/D,, where 


| n/2 n/2 +++ n/2 
ijn/2 n «++ n/f/2 





i Na = 





in/2 n/2 n 


(n — 2) columns 





This determinant may be evaluated by subtracting the first row from 
each of the other ones and expanding the resulting expression by minors 
along the first column and continuing with expansions of the same kind 


P : ‘ ° — . N \n-2 
successively with the resulting expressions, obtaining finally {= }"~. Thus 
m, b 2 J 


Na = -(3), and Cs comes out to be —2/n’. The determinants 


from which Nie, Noa, +++ , ANd N(n_2)(n—-1) are to be evaluated will differ 
in every case from that used to evaluate N.», but by proper inversions of 
rows and columns they can be brought to coincide with the above one, so 
that Noe = Nog = -*+ = Nn-2)~n-1) and therefore, Ca = Cae = -°* = 
C(-2)(n—1) = —2/n’. 

In experiments of this nature where each treatment is repeated 2(n — 1) 


times, Nii is (n = 1) Gyr = (n bas 1)(n)" . Ni; = -(3)2 Es 


-(n)"’, and D, = o(3)@™ = n(n)". Therefore, Cis = 
“a I Par  * as n—2 a ; 
( ae “a ._— , and Ci; = a, = ss That is, as the sums 


of squares and products double, the values of the C;;,’s and C;,’s are divided 
j 
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by 2. If the treatments are repeated 3(n — 1) times, Ci: = (20 = 0) 
n 
and C;; = (=) , and similarly for other multiples of n — 1 replications, 


In the above discussion no reference has been made to C,, and to the 
C'i,’s, that is, the values corresponding to the last variety. Now C,,, = 


GC. = Cw Ss eee = C(n—1)(n—1) and Cas = G5 = Ci. sa eee = C(n—2)(n-1)« 

This may be proved as follows: VN = —(A + B+ C+---+ (N —1)), 

and therefore, for n — 1 replications, Can = —(Caa + Cas + Cac +++ 

+ Cam-1) _ —Cuaa es (n rs 2)Cw = ae > 2) ae (n - 2) (3) = = - 
n n n 

Ca = Cac = +++ = Cin-2)—n-1). This holds also for Con, Con, -°* » Cinna. 
’ Y ~Z 

Also, On = — (Can + Con + gi 2s + Cuxte) a —(n i 1) (+) = 


2(n — 1)/n? = Caz, etc. The above relations hold also for other mul- 
tiples of (n — 1) replications. 


Numerical example 


Diagram II shows the geographical distribution of the plots, together 
with the treatments and yields per plot, in cwts. per acre, of a fertilizer 
experiment with cotton performed by Messrs. A. Riollano and J. Pastor 
Rodriguez in cooperation with the author, on Coto clay at Mantilla Farm, 
Isabela. 

Column (1) of table I shows one way in which the plots may be paired 
so that one may employ the procedure of calculation discussed above. 
As will be noticed, the paired plots lie side by side in all but two cases: 
where plot 56 has been paired to plot 46, and where plot 35 has been paired 
to plot 47. In these two cases the topography of the field in the section 
where these plots were located suggested that they should be paired. 
Columns (2) to (8) show the coefficients of the treatment constants cor- 
responding to each of the paired plots. The constant of treatment H 
has been expressed in terms of the constants of the other seven treatments 
according to the relation H = -(A+B+C+D+E+F+4+60). 
Column (9) shows the respective differences between the yields of the paired 
plots, taken in the order indicated in. column (1). Though unnecessary, 
the plots of each pair were paired so that these differences in yields were 
positive. The numbers in the central portions of columns (10) to (16) 
are the products of the entries in columns (2) to (8) by those of column 
(9), as expressed at the tops of the respective columns. The first of the 2 
rows at the bottom of columns (10) to (16) consists of the sums of the 
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entries in the body of each of the respective columns. 
second row at the bottom of these columns are each one-half of the entries 


of the first row and indicate the values of Say, Sby, --- 





DraGcram II 
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The entries of the 
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Geographical distribution of plots of fertilizer experiment with colton 
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to whether the heading of their respective columns are AY, BY, --- , GY. 
The theoretical reason for this procedure of finding Say, Sby, --- , Sgy 


has been given elsewhere (1). 


Now, in this experiment each of the 8 treatments has been paired once 
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with each of the other treatments, and as mentioned in the section on the 
theory of the procedure of calculation, 


Coa = Cw = °** = Cog = Cu = £2 (8 — 1)/8° = 14/64 = 7/32, 
and 
Ca = Cac = 1+ = Cre = Cra = +++ = Cy = —2/8° = —2/64 = —1/82. 
Knowing, then, the C;;’s and the C;,’s, together with Say, Sby, --- , Sgy, 


one can calculate the yield constants corresponding to each treatment and 
test the significance of differences between these constants. 
Thus, 


A= ba = CaaSay + CapSby + CacScy + CaaSdy 
+ CaeSey + CasSfy + CagSgy 

(7Say — Sby — Scy — Sdy — Sey — Sfy — Sgy)/32, 
(8Say — (Say + Sby + Scy + Sdy + Sey + Sfy + Sgy))/32 

= (—0.8 — 48.00)/32 = —48.8/32 = —1.525. 
B = bs = (8Sby — 48.00)/32 = (70.00 — 48.00)/32 = 22/32 = 0.6875 
C = be = (8Scy — 48.00)/32 = (44.80 — 48.00)/32 
. = —3.2/32 = —0.1000 
D = bp = (8Sdy — 48.00)/32 = (50.80 — 48.00)32 = 2.8/32 = 0.0875 
E = by = (8Sey — 48.00)/32 = (72.40 — 48.00)/32 = 24.4/32 = 0.7625 
F = br = (8Sfy — 48.00)/32 = (76.40 — 48.00)/32 = 28.4/32 = 0.8875 
G = be = (8Sgy — 48.00)/32 = (70.40 — 48.00)/32 = 22.4/32 = 0.7000 
H=-(A+B+C+D+£+4+F+4G) = —1.5000. 


To make comparisons between these statistics, one must calculate their 
standard errors. The sum of the squared deviations of the yields corrected 
for differences in fertility between the different 2-plot blocks is found by 
dividing the sum of the squares of the entries in column (9) of table I by 2. 
The squares of the corresponding entries constitute the body of column (17) 
of table I, their sum is the next to the last item in that column, and the 
required sum of squared deviations corrected for differences in fertility be- 
tween the 2-plot blocks is one-half the next to the last item in column (17), 
that is, 57.18. This sum of squared deviations is subject to 28 degrees of 
freedom, since 28 degrees of freedom have been lost: one df corresponding 
to the mean of all the yields and 27 df due to the fitting of the 27 different 
block constants necessary to express the relations between the fertilities 
of the 28 2-plot blocks. 

The reduction of the sum of squared deviations due to the fitting of the 
treatment constants to the yield data is, according to the usual relation, 
Sy'y’ = ASay + BSby + --- + GSgy = (—1.525)(—0.1) + 0.6875 
(8.75) + --- + 0.7000 (8.80) = 27.70, subject to 7 df. 
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The reduced sum of squared deviations is, therefore, S(Y — Y’) = 
Syy — Sy’y’ = 57.18 — 27.70 = 29.48, subject to 28 — 7 = 21 df. The 
estimate of the reduced variance, V’, is then 29.48/21 = 1.4037. 

The reduction in the sum of squared deviations due to the use of the 7 
treatment constants may be tested by means of the F-test as follows: 
F = 27.70/7(1.4037) = 27.70/9.8259 = 2.82. Since the value of F at the 
5% point, for 7 and 21 df is 2.49, the reduction in the sum of squared de- 
viations due to the fitting of the treatment constants is a significant one, 

The total sum of squared deviations of the yields of this experiment is 
134.29, and as found above, the sum of squared deviations of the yields 
corrected for differences in fertility between blocks is 57.18. The fitting 
of the block constants led therefore to a reduction of 134.29 — 57.18 = 
77.11 in the sum of squared deviations. These relations are indicated in 
table IT. 














TABLE II 
Analysis of the total sum of squared deviations following the 2-plot block method 
a ogee ere eT | 
Sumas | F values 
St a 2) | Degrees ie Variance | : 
Source of the deviations FN yp Pr orl eatinate maak 
+ mental 5% 1% 
MMS. ire oh carck. soe ol 55 134.29 | | 
Blocks 27 ye) | 
Treatments..........:.. 29.48 4.21 2.82 | 2.49 3.65 


27.70 1.4037 | 
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In testing by means of the t-test any difference between the fitted treat- 
ment constants, the variance of any such difference would be found by 
multiplying V’ by the corresponding factor of the nature of Ci; — 2Ci; + 
C;;. Since, in the case of this experiment all the C;;’s are equal and all the 
C;;’s are equal, the variance of the difference between any two treatment 
constants is (7/32 — 2(—1/32) + 7/32) 1.4037 = 0.70185. The standard 
error of any such difference is, therefore, (0.70185)"” = 0.8378, and any 
such difference, to be significant at the 5% point, must exceed 
2.080(0.8378) = 1.74 ewts. cotton per acre. 

The treatments compared in this test are described in table III. The 
corrected mean yield of any given treatment is found by adding algebraic- 
ally the respective treatment yield constant as found above to the mean 
yield of all plots, 6.64. 

Since treatments E, F, and G differed only with respect to the rate of 
the K.O applications, the comparisons of their corrected mean yields or, 
what is equivalent, of their treatment yields constants, will furnish in- 
formation on the effect of the K,O applications on the cotton yields. Such 
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comparisons indicate that the effect of the K,O applications on the cotton 
yields was not significant since the largest difference between any two of 
said constants was only 0.1875, whereas such a difference would have to 
exceed 1.74 ewts./A. to be significant at the 5% point. 

Treatments C’, D, and G varied only with respect to the P.O; applications. 
In this case also, none of the differences between the yield constants were 
significant although they were of a larger magnitude than in the case of the 
K,0 applications. 

A similar comparison of the yield constants of treatments A, B, and G 
indicates that there was a significant increase in yield, caused by the ap- 
plication of at least 100 lbs. NHsz per acre, since the difference between the 
yield constants of treatments A and B was 2.22 cwts. cotton per acre, and 
the corresponding difference between the yield constants of treatments 


TABLE III 











Description of treatments and corrected mean yields 
| Lbs. of nutrient applied per acre | ygeatment Pengo 
| NH: P20s | K:0 | (cwts./acre) 
| 

A eo | =~ | @ | =a 5.11 

B a ee 7.33 

C 200 | 0 200 —0.1000 6.54 

D 200 | 100 200 0.0875 6.73 

E 200 | 200 0 0.7625 7.40 

F 200 | 200 100 0.8875 7.53 

G 20 | 20 $| 200 0.7000 7.34 

H  S 0 | 0 —1.5000 5.14 











A and G@ was 2.23 ewts. cotton per acre, both differences being significant 
at the 5% point and nearly so at the 1% point. The difference between 
the yield constants of treatments B and G was not significant indicating 
that the second 100 lbs. NH; applied per acre did not affect significantly 
the yields already obtained With the first 100 lbs. NHs applied per acre. 

The yield constant of treatment H was almost equal to that of treatment 
A, corroborating the conclusions drawn above that nitrogen was the only 
one of the three elements tested which affected the yields significantly. 

This experiment might have also been interpreted by using only 3 con- 
stants to explain the variations in yield caused by the various treatments 
as described in the article (1) previously referred to. Such a study was 
made yielding the same results previously obtained, except that the partial 
regression coefficient of the yields on the nitrogen applications came out 
to be significant at the 1% point. 
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SUMMARY 


A new method of performing field experiments with relatively sms 1 
numbers of treatments is described. The requirement to be fulfilled by the® 
layouts of such field tests is specified and examples of possible designs for a 
5-treatment experiment are illustrated. The theory of the procedure of. 


calculation is discussed 


furnished in connection 
performed with cotton. 


and a numerical example of said calculations ig) 
with the interpretation of a fertilizer experiment: 
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